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Abstract 

We propose an alternative to the usual time-independent Born-Oppenheimer ap- 
proximation that is specifically designed to describe molecules with non-symmetrical 
hydrogen bonds. In our approach, the masses of the hydrogen nuclei are scaled dif- 
ferently from those of the heavier nuclei, and we employ a specialized form for the 
electron energy level surface. As a result, the different vibrational modes appear at 
different orders of approximation. 

Although we develop a general theory, our analysis is motivated by an examination 
of the F H Cl~ ion. We describe our results for it in detail. 

We prove the existence of quasimodes and quasienergies for the nuclear vibrational 
and rotational motion to arbitrary order in the Born-Oppenheimer parameter e. When 
the electronic motion is also included, we provide simple formulas for the quasienergies 
up to order that compare well with experiment and numerical results. 



*Partially Supported by National Science Foundation Grant DMS-0600944. 
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1 Introduction 



This is the second in a series of articles devoted to the study of vibrational levels associated 
with hydrogen bonds. The first paper [S] deals with stretching vibrations of the hydrogen 
bond in the symmetric case in which the hydrogen binds two identical atoms or molecules. 
Our prototypical example is FHF", which displays strong anharmonic effects, coupling 
between vibrational modes, and a low frequency for the vibration of the hydrogen along 
the F-F axis. This second paper deals with all the vibrations and rotations in the non- 
symmetric situation. Our canonical example is F H Cl~, which displays weaker anharmonic 
effects and a high frequency for the vibration of the hydrogen along the F-Cl axis. 

Both of our papers contain two main new ideas. The first is the same for both papers. 
Standard Born-Oppenheimer approximations keep the electron masses fixed while all the 
nuclear masses are taken proportional to e~^. We take the hydrogen mass proportional to e~'^ 
while keeping the heavier atoms' masses proportional to e~^. This is physically appropriate 
for many molecules of interest: If the mass of an electron is 1 and e is defined so the mass 
of a carbon C^^ nucleus is e~^, then e = 0.0821, and the mass of a nucleus is 1.015 e~^. 

The second novel idea is to exploit the smallness of certain derivatives of the electron 
energy level surface for the molecule being studied. Here our two papers are completely 
different, and they are motivated by examinations of numerically computed electron energy 
level surfaces using Gaussian 2003 software [3]. In the symmetric case, the second derivative 
associated with moving the H along the axis oi AH Ais small, and we could allow it to be 
small and negative if the H nucleus felt a double well potential. In the non-symmetric case, 
if the H is more weakly bound to the B m A H B , we assume all the derivatives associated 
with moving the B relative to AH in AH B are small. We assume all derivatives associated 
with stretching the distance between A and H not to be small. 

To describe the smallness of the small derivatives, we could have introduced another 
small parameter. Instead, we have elected to let e play a second role. We take all the small 
derivatives to be proportional to e. For the choice of e = 0.0821 indicated above, that is 
again appropriate for our F H F~ and F H Cl~ examples. The small derivatives are on the 
order of e in units where the non-small derivatives are on the order of 1. 

We shall now restrict our attention to triatomic non-symmetrical hydrogen bonded 
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molecules AH B, and assume the H is more strongly bound to the A. We do an asymptotic 
expansion for small e, and our main results are the following: 

1. To their respective leading orders, the vibrational levels are described by three inde- 
pendent harmonic oscillators in appropriate Jacobi coordinates: two separate one- 
dimensional harmonic oscillators and one two-dimensional isotropic harmonic oscilla- 
tor. This is in contrast to the usual Born-Oppenheimer theory in which one obtains 
one coupled four-dimensional harmonic oscillator. Our technique does not require go- 
ing through the diagonalization process to separate the normal modes. The different 
modes appear at different orders of the expansion, in contrast to the Born-Oppenheimer 
situation, where all vibrations are of order e^. 

2. The highest frequency vibrational states have energy of order e^/^. These are the stretch- 
ing oscillations of the A-H bond with the B approximately sitting still. 

3. The next highest frequency vibrations are the two degenerate bending modes. They are 
of order e^. 

4. The lowest vibrational energies are of order e^/^. They are the stretching oscillations of 
the weak bond between the A H and the B. 

For the specific case of FHCl^, we have the following comparison of results, where 
vibrational energies are measured in cm~^. The experimental results come from [2]. We 
note that the experiments were not done in the "gas phase," so they may not accurately 
represent results for the isolated ions. All the Gaussian 2003 results presented in this paper 
are obtained by using the MP2 technique with the aug-cc-pvdz basis set. The software 
implements the standard Born-Oppenheimer approximation. The results for our model 
come from approximating the ground state electron energy surface with Gaussian 2003 and 
then applying our techniques. 



Mode 



Experiment 



Gaussian '03 



Our Model 



F -H stretch 



2710 



2960 



2960 



bends (degenerate) 



843 



875 



871 



FH - CI stretch 



275 



246 



251 
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Remarks 

1. It is not surprising that the results for our model are close to those obtained by Gaus- 
sian since we have used the same electron energy surface. The Gaussian software deals 
with the full 4-dimensional harmonic oscillator, whereas our technique deals with two 
1-dimensional harmonic oscillators and one isotropic 2-dimensional harmonic oscilla- 
tor. Evidently the Jacobi coordinates we have chosen are very close to the normal mode 
coordinates for the 4-dimensional oscillator. 

2. The results from Gaussian and our model are just leading order (harmonic) calcula- 
tions. Including higher order terms from the expansions might bring these into better 
agreement with experiment. Also, we again emphasize that the experimental results 
were not obtained for isolated ions. 

A recent chemistry article [9] contains data for vibrations of eighteen hydrogen bonded 
molecules in the gas phase. It also contains an idea for quantifying how symmetric or non- 
symmetric a hydrogen bond is. Its conclusions are consistent with the analysis in our two 
papers. Figure 2 of that article plots the vibrational frequency of the A — H stretch versus 
the difference in the "proton affinities" of A and B for a molecule AH B. When A and 
B are identical, the frequency is low (800-1000 cm~^), and when they attract the proton 
very differently, the frequency is high (1600-3500 cm~^). In our symmetric analysis, this 
vibrational energy is of order e^, whereas in our non-symmetric analysis, it is of order e'^/^, 
which is roughly 3.5 times larger when e = 0.0821. 

Remarks 

1. We assume that the ground state electron energy level we are considering is non- 
degenerate for all nuclear configurations of interest. Thus, we do not consider situations 
that exhibit the Renner-Teller effect [HI [101 E] • 

2. Since our analysis includes rotations of the whole molecule, some small effects show up 
in the calculations. For example, /-type doubling [7] occurs for terms that have non- 
zero eigenvalues of the L^i operator at low order. [L^i is the nuclear angular momentum 
around the A—B axis.) States corresponding to L^' eigenvalue ±k with k > 1 generically 
have their degeneracy in energy split at order g^+sfc q^j, model. 
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The paper is organized as follows: In Section [21 we describe our model in detail. In 
Section [3l we do the semiclassical expansion to all orders for the nuclei. In Section |4] we 
include the electrons. However, when we include the electrons, we just show that the energy 
expansion is valid through order e^. Going to higher order is extremely complicated. 

2 Semiclassical Analysis for the Effective Nuclear 
Hamiltonian 

In this section, we give a precise description of the Hamiltonian for the nuclei. As mentioned 
above, we consider a molecular system AH B in which the hydrogen is much more tightly 
bound to the A than to the B. 

We construct the coordinate system we use in two steps, as illustrated in the figures 
below. The first step is to choose a standard Jacobi coordinate system for the nuclei in their 
center of mass frame of reference. The first three coordinates are the components Xi, X2, 
and X3 of the vector X from the A nucleus to the H nucleus. The fourth, fifth, and sixth 
coordinates Yi, Y2, and Y3 are the components of the vector Y from the center of mass of 
the A and H nuclei to the B nucleus. 



z 




Jacobi coordinates for the molecule 
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We now change from these coordinates to new ones that we call {Y, 6, (f), R, 7, X). The 
(Y, 6, 0) are spherical coordinates for the vector described by {Yi, Y2, Y^) in the original 
center of mass frame of reference. The {R, 7, X) are cylindrical coordinates for the vector 
{Xi, X2, X3) in a frame of reference that rotates so that the axis for these coordinates is in 
the direction of the vector described by (Yi, Y2, la). The precise definition is below. 




Jacobi coordinates fixed at the origin 
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The final coordinate system 



One reason for using tliese coordinates is tliat tlie potential energy surface depends only 
on Y, X, and R. A second reason is that in these coordinates, we can separate the total 
angular momentum and its z component from the other motions easily. Also, to low 
order in perturbation theory, the angular momentum L^' conjugate to 7 (which is the angular 
momentum in the direction of {Yi, Y2, Y3), gives another convenient quantum number. Note 
that Lz' does not commute with the full Hamiltonian. 

The drawback to using this coordinate system is that the kinetic energy expression is 
quite messy. The complication comes from the Laplacian in the (Y, 6, 0) variables. The 
Laplacian in [R, 7, X) is simply the usual cylindrical Laplacian. 

These coordinates are closely related to ones used in [1] to deal with Born-Oppenheimer 
approximations for diatomic Coulomb systems. There is a minus sign error in the expression 
for L ■ J term on page 32 of that paper. 

As mentioned above, (Y, 6, 0) are just standard spherical coordinates. To describe the 
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other three coordinates precisely, we first define the rotation 

/ cos(6') cos(</)) — sin(0) sin(6') cos(0) \ 
cos(6') sin(0) cos(0) sin(^) sin(0) 
- sm{e) cos{9) 



It maps the vector 



coordinates (G, 6, 6) by 




V 



to the unit vector in the direction of 




We then define 




7^l( 



1-1 




Next, we define another rotation 



7^2(7) 



/ cos(7) — sin(7) \ 
sin(7) cos(7) , 
1/ 

where, for generic vectors C,, 7 is defined by requiring the second component of [7^2(7) C, 
to be and its first component to be positive. We then define coordinates X and R by 



:^2(7): 





Our Hamiltonian has kinetic energy 

^3 



A 



(Xi, X2, x^) 



A 



iYi,Y2,Y3), 



2/ii(e) — — 2/i2(e) 

where fii{e) and /i2(e) are modified reduced masses that we describe in detail below. Since 
Laplacians are rotationally invariant, under our coordinate changes, the first term simply 
becomes the usual cylindrical Laplacian 



g2 



1 A 

R OR 



1 



2/ii(e) \dR^ R dR R^ d-f^ dX^ 

By a very tedious calculation, the second term in the kinetic energy is 

^ 2d 
^ Y W 



where 



J' 



2/i2(e) 
86^ 



- ^ {J'-2L.J + L^}^ 



1 



sin' e 



92 



92 



dr 



2 cos e 



sin' e 



(2.1] 



is the total angular momentum operator, 

T 1 ( 1, ■ ^ If ■ ^ ^ d \ ( I d d 
/ d ^ d X . d\ d 

+ ^^^°^7^-^cos7^ + ^sm7^J ^ 

and 



dXdR dR^ i?2 V R J dR dX d-/^' 



The modified reduced masses are fiAe) = — ^ and 

/i2(ej = e — , where the three nuclei have masses e m^, e ttlh, 

and e'^^niB- The modified reduced masses have limits as e tends to zero. To isolate the leading 
behavior, we abuse notation and define ui = lim uife) = mu and uo = lim uoU) = 

ihen we have 



Similarly, 



e3 



2/ii(e) 2/ii 2 mA 



^4 ^4 ^5 



2/i2(e) 2/i2 2mA (m^ + 2em//) ' 

We define the operator 

g4 ^5 

= " 2;;^ + 2mA(mA + 26mH) 

written in the new variables, so that the kinetic energy can be expressed as 

^3 g4 

all written in terms of (Y, 6, 0, R, 7, X). 

The quantum fiuctuations of the nuclei around their equilibrium positions occur on short 
length scales, so we now do the appropriate rescaling of variables. We assume the ground 
state electron energy surface has a minimum at F = Yq, R = (because the Hydrogen bond 
is linear), and X = Xq. Under the rescaling, the angles 6, (p and 7 remain unchanged, but 
we replace Y, R, and X by 

y = (Y-Yo)/e^/\ r = R/e^/\ and x = {X-X^)/€''\ 
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Under this rescaling, the total kinetic energy operator becomes 

3/2 Q2 f Id 1 92 \ 6^/2 92 , , 

+ Zir + -3.1^.] - TT-TTTo (2.2) 



2 /xi c?a;2 2 /ii \ dr"^ r dr d'y'^ ) 1^2 dy"^ 

gl3/4 Q ^4 

" ^2 {Yo + eV^y) 9^ + 2/., (Fo + e^'^yY [j'-^L-J + L^] + Die), 
where is still given by (12. ip . but L-J and are now given by the e-dependent expressions 



L-J = ( r sin 7 I- - e~'/'iXo + e'/^x) sin 7 ^ - e^^/^ + ^ 
\ ox or r 



, I d ^ d 

X I 7:— — cot ^ — 

sin 6 d(j) 07 



+ ( e-^/V cos 7 I- - e-"\X, + e'"x) cos 7 |- + e-V2^^±i!^ ^ 
\ ax ar r 07 / 
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and 



L' = - e-"'r' 1^ + 2 e-'" (Xo + r - e-^(Xo + e^^/^x)^ 

ax^ ax ar ar"^ 

+ e er ' 



(97^ Y r J dr 



+ 2e-3/4(Xo + e3/4x) ^ 



(?x 972 



Remarks 

1. The operator L ■ J can be rewritten as 



-1/4 L ^ (t - t \ - .-1/2 ^+e^/^a; 9 



where 

ab' 07 sm 6 04. 
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By explicit computation, one can verify that L+/ and L„/ commute with both and 
Jz- The operators L+z and L_i are raising and lowering operators for the eigenstates of 
Lz'- 

2. The dominant order terms in the expressions in L ■ J and are the ones of order in 
. Because of the overall factor of that multiplies these operators in the Hamiltonian, 
they are not relevant until the order perturbation calculations. 

Motivated by numerical calculations for the FHCl~ ion, we assume the ground state 
electron energy surface near its minimum depends only weakly on R and Y . To exploit this, 
we decompose the potential energy surface as 

V^{X) + eV2{X,R,Y), (2.4) 

where Vi and V2 have Taylor expansions of the forms 

00 

Vi{X) ~ ao + ^ aj^X-XoY, and (2.5) 
i=2 

V2{X,R,Y) ~ b,,k,i{X-XoyR'{Y-Yoy. (2.6) 

j + k + l>2 
k + l>l 
k even 

The restrictions on the indices in V2 are obtained requiring all pure X dependence to be Vi 
and by requiring V2 to be even in R (because of the symmetry). 

We now can state our results for the semiclassical analysis of the bound states for the 
nuclei. 

Theorem 2.1 Consider the Hamiltonian 

Hie) = - A(x„x„X3) - ^iYuY„Y,) + V^{X) + e V2iX, R, Y), 

rewritten in terms of the variables {X, R, Y, 6, (p, 7). Assume Vi and V2 are functions 
that satisfy Ii2.5\) and Ii2.6\) . Assume V\ has a unique global minimum at X = Xq > 0, 
with a2 > in i\2.5\) . and that liminf|x|^oo V{X) > ao- Assume V2 has a unique global 
minimum of at X = Xo, R = 0, and Y = Yq > 0, with 60,2,0 > and 60,0,2 > in 

N 

112. 6\) . Given any integer N > 0, there exist a quasimode "^w/iie) = e^^^ipi/^ and a 

1=0 

11 



N 



quasienergy E]\f/4{e) = ''^^ ^''^'^ £i/i, such that Wi^i/il 

1=0 

each I, and 

I (^H{e) - ^7v/4(e)) ^N^e) | 

for some Cn that depends on the choices of n, k, m, and p below 
Furthermore, 



= 0(1) for each I, Si/^ = 0{1) for 



^0 — flO; ^1/4 — ^2/4 — ^3/4 — ^4/4 — ^5/4 — ^7/4 — ^9/4 — ^11/4 



6/4 



a/ 2 02//^! ( + \ 



for n = 0, 1, 



^8 



2feo,2,o//^i (2m + \k\ + 1), for an integer k, and m = 0, 1, 



10/4 



2feo,o,2//"2 { P + 2 



for p = 0, 1, 



and £i2/4 is given by the expression ^3. 7| j. The rotational energy first appears in Eie/^. 
For fixed angular momentum quantum numbers j and j^, for order N > 12, the states with 
k = are non- degenerate, and the states with \k\ > have multiplicity at most 2. 



Remark Theorem 12.11 is stated with global hypotheses and without growth conditions on 
the potential. When the electronic motion is also included, the potential energy surface may 
only exist locally. The cutoff functions that are introduced in Proposition 13.21 allow us to 
obtain analogous results with only local assumptions. 

For the FHCl" ion, we have calculated values for the first few coefficients in the expansion 
for V, based on numerically differentiating results from Gaussian 2003. Here distances are 
measured in Angstroms, energies are in Hartrees, and we have used e = 0.0821. 

ao = -560.160 

a2 = 0.567 
^0,2,0 = 0.597 
6^0^^ = 0.853 
bo,o,2 = 0.664 

The e in 02.41) refiects the weakness of the hydrogen bond, and also that the molecule can 
bend easily. The FHCl~ ion essentially looks like a slightly deformed FH molecule with a 
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CI ion quite a long way from the FH. Gaussian 2003 assigns charges associated with each 
atom, and it obtains: 

F -0.58 
H 0.51 
CI -0.93 

The calculated F-H distance is 0.98 Angstrom, and the H-Cl distance is 1.91 Angstroms. 
(For HF alone, the charges are ±0.33, the H-F distance is 0.925 Angstrom, and the calcu- 
lated vibrational frequency is 4083 cm~^.) 

Experimental values [2] for the vibrational frequencies of FHCl~ (in cm~^) are 

275 FH oscillates relative to the CI 
843 bends (2 degenerate modes) 
2710 FH oscillates 

Gaussian 2003 calculates the harmonic vibrational frequencies (in cm~^) to be 

246 FH oscillates relative to the CI 
875 bends (2 degenerate modes) 
2960 FH oscillates 

To leading order, our model has these frequencies proportional to e^/^, e^, and e^^^ respec- 
tively. The specific harmonic frequencies that we obtain for FHCl~ are 

251 FH moves relative to the CI 
871 bends (2 degenerate modes) 
2960 FH oscillates 

3 The Perturbation Expansion for the Nuclei 

We now do the perturbation expansion for the semiclassical motion of the nuclei under 
the global hypotheses of Theorem 12.11 When the hypotheses are satisfied only locally, see 
Proposition 13. 2[ 

The perturbation expansion describes the small e dependence of the eigenvalue problem 
for the following differential operator 

^3/2 Q2 f Id I \ 6^/2 ^ ^ 

+ zir + -.1^.] - (3-1) 



2 /ii dx^ 2 111 \ dr"^ r dr / 2/^2 dy"^ 

,13/4 g A 



' ' { J2 - 2L- J + L^}, 



/i2 (Fo + e3/4y) dy 2^l^{Yo + e^l^yf 



i=2 j ^k + l>2 

k + l>l 
k even 
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At this point we should make the Ansatz that the eigenvalue and eigenfunction have 
expansions of the forms 

oo oo 

£ = ^ i^,(e) Sq^ and ij{x, r, y, 9, 0, 7) = ^ ^^{e) ^pq.ix, r, y, 9, 0, 7). 
/=o 1=0 

Here, 1^0 (f) = 1, "^o is non-trivial, and ui^i^e) / ^[{e) — > as e ^ 0. 

However, one learns that every z/;(e) that occurs is some power of e^^^, so it is somewhat 

simpler just to take i^i(e) = e'/"^, i.e., 

00 00 
£ = ^ e'/^ Si/^ and i/j{x, r, y, 9, 0, 7) = ^ e^/^ ipi/4{x, r, ?/, 9, 0, 7). 
;=o z=o 

Our Hamiltonian, J^, and all commute with one another, so we can simultaneously 
diagonalize these three operators. The eigenvalues of are j(j + 1), where j = 0, 1, 2, . . ., 
and for a given j, they have degeneracy (2j + 1)^. We henceforth use the specific basis for 
the eigenspace for fixed j that is given in Section 4.7 of pQ: 

{ I 3, 3z, k) : = -j, -j + 1, . . . , j; k = -j, -j + 1, . . . , j }, 

where 

Jz I j, j^, /c) = jz I j, j^, /c) and | j, j^, /c) = k \j, j^,, k), 

d d 

where Jz = —i tt- and Lz> = —i — . Note that although J^, J^, and Lz' all commute with 
(70 a7 

one another, L^/ does not commute with the Hamiltonian. 

For future reference, we note also that the operators in (12.31) have 

L+' \j, jz, k) = a+jj^,fc \j, jz, k+1) and L_/ \j, jz, k) = \j, jz, k-1), 

for some When \k\=j, a^jj^j = and _j = 0. 

By restricting attention to given values of j and jz, the wave functions in our expansion 
can now be regarded (with some abuse of notation) as 

j 

i^i/iix, r, y, 9, 0, 7) = ^ i^i/iix, r, y, k) \j, jz, k). 

k=-j 

We now substitute the Ansatz into the eigenvalue equation and equate terms order by 
order. We do not worry about normalization, but produce a quasimode that is 0(1) as e 
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tends to 0. To simplify some of the discussion, we take orthogonal to ipQ for / > 0. The 
results of these computations yield the formal expansions of Theorem I2.1[ 



Order e° These terms simply require oq ^Jq = Sq ipQ. So, 



Order e'/^ for 1 < 1 < 5 The terms of these orders successively require 

Oo ^Z/4 = ^0 -ipl/A + ^1/4 -ipO- So, 

= 0. 

Order e^^^ These terms require — — ^ + 02 i/jq = E^u ipQ. 



This forces 



and 



2/ii dx"^ 

= + 2^ a/ 2 02//^! for some n = 0, 1 
?/;o(x, r, /;;) = /o(r, y, /c) 



where 



$i(x) = (2a2/ii)'/^ vr"i/^ 2-"/^ (^!)-i/2 ^^(3,') g-'Vs 
with x' = (2 02 /ii)^^^ x. The function /o is not yet determined. 



Order e^^^ We introduce the notation 



1 92 



0,1 



+ a2 X 



2jii dx"^ 

Then the e'^^^ terms require [-ffo.x — £e/i\ "^1/4 = ^7/4 V'o- 

We first examine the components of this equation that are multiples of These 
components require 

^7/4 = 0. 

We then examine the components that are perpendicular to $i(x) in the x variables. These 
±a; components require 

i^i/4{x, r, y, k) = /i/4(r, y, k) 
15 



where the function /1/4 is not yet determined. 



Order e^^^ These terms 



require 



[i^o,x - ^6/4] V'2/4 - ^(^^ + -^ + ^^J + ^0.2,0 r ^0 = ^8/4 ^0- 



The 1 1 a; components of this equation require 

-f^0,r,7 V'O = ^8/4 V'O, 

where 



1 , 1 9 , 1 92 ^ , , 2 



0''"'')' 2jJ,i \ dr"^ ~^ J, Q^. ^ j ~^ 

This is a standard isotropic two dimensional Harmonic oscillator problem that one can solve 

dip 

by separating variables. In our context, the angular operator L^/ = — i^ — has eigenvalues 
k = 0, ±1, ±2, ■ ■ ■ , ±j and eigenfunctions e*'^'^. For each such k, the radial operator 

2/ii y (?r2 r 9r y 

has eigenvalues 

^8/4 = {2 m + \k\ + 1) ^^2 5o,2,o/Ati, where m = 0, 1, ■ ■ ■ . 
The corresponding normalized eigenfunctions are 



/ ^ ("^0 /o/, „ n1/4 / n|fc| r|fc|/ /2s ,.'2/2 

V(m+|fc|)! (2&0'2,o^i) (r) 4n (r )e , 

where, r' = (2 69,2,0 /^i)^''^ '"j m > 0, and Lm is a Laguerre polynomial. 

We permanently fix one such value of Ss/4. Since different pairs (m, k) can occur, we 
define 

K = {keZ : \k\<j, and m{k) > }. 

where 

1 



m{k) = - l^Ss/i j y2 6o,2,o//ii - \k\ - \ 

One can easily show that K is non-empty and has at most j + 1 elements. 
For k G we define the normalized wave functions 



^2{\k\, r) 
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and take 

/o(r, y^k) = < 

[ otherwise. 

The functions go{y, k) ioi k & K are not yet determined. However, we now have 
ipoi^x, r, y, 9, 0, 7) = ^ go{y^ k) <l>i(x) ^2i\k\, r) | j, j^, k) . 

For future reference, we let Zi denote the subspace spanned by 

{^i{x) ^2{\kl r)\j,j„k) : keK}. 

The -Lx terms at this order require [-f^o,x — ^e/i] ^2/4 = 0, which simply forces 

^2/4 = /2/4(?^, y, k) $i(a;). 

Order e^/^ These terms require 

[i?o,x - i'3/4 + [Ho,r,y - £8/4] + as ^0 

= £^9/4 i^o- (3.2) 
The 1 1 a; components of this equation require 

[Ho,r,^ — £8/4] ^1/4 = £^9/4 V'O- (3.3) 

We first examine the components of this equation that belong to the subspace Zi. These 



X Wz-^ components require 



£^9/4 = 0. 



fi/4{r, y, k) 



Next, the J-Zi components of (13. 3p that are orthogonal to Zi require [i?o,r,7 — £^8/4] ^1/4 = 
0. This forces us to choose 

9i/4{y, k) $2(|A;|, r) li k e K 
otherwise. 

The ^-x components of ( 13. 2p require [Hq^x — £^6/4] ips/i + o-s i^o = 0. We solve this 
equation by applying the reduced resolvent operator [Hq^ — £q/a\~^ ■ The result is 

^3/4(2:, r, y,k) = - as fl'o(z/, k) ^2i\k\, r) [Hq,^ - £&/A]r^ {x^ ^i{x) ) 

+ /3/4(r, y, k) $i(x). (3.4) 
17 



Order e^"/^ 



[Ho^x — ^6/4] i^A/i + [-f^0,r,7 — ^^8/4] V^2/4 



1 d'^lpQ 

2/7 "9^2" + '^3 x'^i^i/i + Vo,2 y'^ ^0 + ^1,0,1 3;y^o = £w/4 V'o- (3.5) 

7^ + 00,0,2 y Wo = ^10/4 

2/i2 (9?/2 



The IIj; \\zi components require — + 60,0,2 Z/ V'o = ^10/4 "^o- This forces us 

to choose 

^10/4 = (^^+0 y^2 6o,o,2//U2 where p = 0, 1, 

and 

^7o(l/, fc) = Co,, $3(1/) if e i^, (3.6) 

where 

$3(2/) = (2 6o,o,2^2)^/«vr-V4 2-p/2(p!)-i/2^^(^')e-yV2 

with = (2 60,0,2 yU2)^/^2/. 

So far, the Co,, in (13. 6p are arbitrary for k G Ji', but we henceforth assume they satisfy 
the normahzation condition 

^ I Co,fc 1^ = 1. 

fcG/< 

For future reference, we let Z2 denote the subspace spanned by 

{$i(x) $2(|fc|,r) $3(2/) \j,J.,k) : keK}. 
The II a; -Lzi components require 

g2/4{y, k) $2(|fc|, r) if /c e 



/2/4('^, y,k) 



otherwise. 



The components require [i/o,x —^^6/4] "^4/4 + 03 x^'?/'i/4 + 61,0,1 xy ipo = 0. 
We apply the reduced resolvent of ifo,x to obtain 

V'4/4(a;, r, y. A;) = - 03 5-1/4(2/, A;) ^2{\k\, r) [iJo,x - i^6/4]r"^ {x^ ^i{x)) 

- 61,0,1 co,fc ?/ $3(2/) $2(|A;|, r) [i7o,x - ^6/4]^^ (a; $i(x)) 
+ /4/4(r, 2/, fc) $i(x). 
Note that the first two terms are zero ii k ^ K. 



Remarks 

1. At this point, we have completely determined ipo, except for the values of Co,fc for k E K. 
Restoring the angular dependence in the notation, we have 

^0 = J] Co,fc $i(x) $2(|A;|, r) ^^{y) | j, j„ k). 

k&K 

Since j and jz are fixed, this is a linear combination of at most j + 1 linearly independent 
states. 

2. As we shall see, the degeneracy generically partially splits at order e^^^^. At that point, 
states with different values of have different energy, but two states with k = ±A for 
A > have the same £^12/4- In terms of the energy, the degeneracy of these two states 
generically splits completely at order e^"''^'^. When A = 1, this splitting has long been 
observed in the spectra of linear polyatomic molecules. It is called /-type doubling [7] . 

3. We have determined the dominant terms for the eigenvalue: 

£0 + e'/' + V^W^ + (2m(A;) + \k\ + 1) ^^2 60,2,0^1 

This quantity does not depend on the quantum numbers j, jz, ot k E K. 
The dominant contribution to the energy from the total angular momentum is 

—75 — , so it enters at order 16/4. 

2/i2^0 

4. Below we impose the condition that every 1^1/4 with / > be orthogonal to the subspace 
Z2. 

5. At the next order, the pattern emerges for how to do all higher order formal perturbation 
calculations. For / > 11, we have the following: 

• the IIj. II^j IIj^ terms determine S1/4, 

• the 11^. \\zi -i-y terms determine the y-dependence of g{i-iQ)/i{y, k) 

• the \\x J-Zi terms determine the r and k dependence of f{i-?,)/A{r, y, k), and 

• the ±2; terms determine the x-dependence of ip{i-6)/4:{x, r, y, k). 

Since the general pattern occurs at the next order, we present full calculations for only 
one more order explicitly. 
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Order e^^/^ 

[Ho,x — ^6/4] V'5/4 + [Ho,r,^ — ^8/4] V'3/4 + [Ho,y — SlQ/4\ 1pl/4 

+ as a;'^^/'2/4 + &i,o,i ^1/4 + ^2,1 y i^o + &i,2,o a;?^^ ^Ao 
The 1 1 a; \\zi \\y terms require 

^11/4 = 0. 

The 1 1 a; \\zi -i-y terms require 

9i/i{y,k) = - 60,2,1 co,fc ($2(|/i;|, r), r2$2(|A:|, r))^, [Ho^y - [y <l>3{y)) . 

for k E K. This is the first place in the perturbation calculations where different values of 
\k\ yield different results. Note that we could add Ci/4^fc ^z{y) to gi/4,{y, k) when k G K, but 
we have chosen ci/4_fc = to impose the condition that ipi/^ be orthogonal to the subspace 
Z2. See Remark 4 above. 
The 1 1 a; -L^i terms require 

[ifo,r,7 - £s/4] h/A + P±z^ Wo.y " £"10/4] /l/4 + &0,2,1 ^ /q = 0, 

where P±z^ denotes the projection onto functions orthogonal to the subspace Zi. We have 
already seen that the non-zero /i/4(r, y, k) belong to the subspace Zi, so P±z^ Wo,y ~ 
^10/4] /1/4 = 0. Thus, applying the reduced resolvent of -^0,^7 (which is zero on Zi), we 
obtain 

h/i{r, y,k) = - feo,2,i Co,fc y $3(2/) [^o,r(|fc|) - ^8/4]^^ *2(|fc|, r) 

+ 93/4{y, k) $2(|A;|, r) li k e K, and 

/3/4(r, y,k) = ii kiK. 
Here, we have used the notation 

i_d_ e_ 

and the direct sum decomposition 

[i^0,r,7 - ^8/4];' = [i/o,r(|A;|)-^8/4];' 

|fc|<i 
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which results from i/o,r,7 commuting with L^/. 
The A-x terms require 

+ as V'2/4 + ^1,0,1 xy ^1/4 + 61,2,0 xr'^ ^jq = 0, 

where i^y^ denotes the component of 'ip-i/4 orthogonal to $i(x) in the x variables. By com- 
bining this with (13 ■4p and (13.61) . we have 

( - Co,k "^siy) '^2{\k\, r) [Hq^x- Se/^]-^ {x^<^i{x)) if keK 
^3/4 ^' y^k) = < 

' y if fc^i^. 

So, we see that [ifo,r — "^3/4 = 0. Thus, we have 

i'5/4ix, r, y, k) = - a-i g2/4iy, k) ^2i\k\, r) {[Ho,x - Se/Alr^ (x^<l>i(x))) 

- bifi^iy gi/^iy, k) <l>2i\k\,r) {[Hq^x - Sq/^;^ 

- 61,2,0 Co,k '^3{y) '^2{\k\, r) ( [Hq^x - ^a/iV '^li^)) ) 

+ h/^ir, y, k) if keK. 

For k^K, 

^5/4(3;, r, y, k) = h/iir, y, k) $i(x). 

Note that only g2/4:{y, k) (for k G K) and frj/i{r, y, k) in these expressions have not yet been 
determined. 

Remarks 

1. Amazingly, 7/^1/4 ^ 0. This component of the wave function involves an anharmonic 
correction related to the bending and AH-B stretching modes. Restoring the angular 
dependence to the notation, we have 

■01/4(0;, r, y, 9, 0, 7) 

= - 60,2,1 Yl ""^'^ (*2(|fc|,r),r2$2(|fc|,r))^ 

k&K 

X ^i{x) $2(|A;|, r) [ifo,s/ - ^10/4];^ {y ^-iiy)) I J, J., k). 

2. Although we do not present the full calculations at order e^^/^, we do calculate £^12/4 
explicitly. It is generically contains non-zero anharmonic corrections. 
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Before going further with the expansion, we present a summary of what has been deter- 
mined so far. 



E = £o + e3/2 \ n + ]^ v^W^ + {2m{k) + \k\ + I) 6o,2,o//Ui 



The last information for £ came from order 11/4 



' I U' 11^1 1 1 J/- 



^0 = XI ^O''^^ *2(|fc|, r) $3(y) \jjz,k). 



k&K 



This was completely determined at order 10/4, H^; ||zi- 

^1/4 = -60,2,1 J] Co,fc {^2{\klr), r2$2(|fc|, r))^ 

X <l>2(r) [ifo,j; - ^10/4];^ (y'^'sly)) |j,J^,A;) 

This was completely determined at order 11/4, ||^-^ 



^^2/4 = 92/i{y: k) ^i{x) ^2{\k\, v) \j,jz,k). 

The last information came from order 10/4, ||^. ±2^. 



^3/4 = - &0,2,1 5Z ^0''^ '^l(^) '^'3(1/)) 

X [HoA\k\) - £s/4];' (Px,^ $2(|A;|, r)) 



The last information came from order 11/4, ||^ J^Zi 
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^4/4 = bo,2,i [Ho^y - Sio/ilr^ (?/$3(y)) [Ho,x - ^e/i]^^ {x^^i{x)) 

Yl co.fe {'^2m,r), r'M\k\,r))^ $2(|A;|, r)| j, j„ ) 



X 

k£K 



61,0,1 (^$3(2/)) [^0,^-^6/4]/ {x $i(a;)) ^ co,fc $2(|A;|, r) I J, j;,/c) 



+ ^ /4/4(r, y, k) $i(a;) | j,j^,A;). 

The last information came from order 10/4, ±3, 
(coupled with 11/4, ||^ Wz^ J-y, because of gi/4). 



V'5/4 = - as Y g2/i{y, k) <l>2{\k\, r) ( [i^o,x - ^6/4]/ (x=^$i(x))) 

- &w W E Co,fc ($2(|A;|,r), r2$2(|fc|,r))^ ( [i/o,. - ^e/J^M^ '^'i(^))) 

X $2(|A:|,r) (2/[i/o,,-^^io/4].T'(y$3(y)))|j,J.,fc) 

- &i,2,o J] co.fc $3(y) r^$2(|A;|, r) ( [ffo.x - ^e/J^^ (x$i(x))) 



+ X] h/4{r, y, k) $i(x) I 
The last information came from order 11/4, ±3;. 



We now return to describing higher orders of the perturbation expansion. We determine 
£^12/4, and explicitly write the equations that must be solved through order e^^/^. That is the 
order at which the angular momentum quantum number j appears, and the degeneracy due 
to rotations is split. 
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Order e^^/^ 

[Hq^x — ^6/4] V^6/4 + [-f^0,r,7 — ^4/4 + [Ho,y — ^lO/i] + ^3 V^3/4 

+ 04 V^o + bifi^ixyip2/i + bo^2,i r'^ y + ^1,2,0 a^^^^ V'i/4 

+ &0,4,0 ^0 + ^^^y2 ^0,r,7 ^/'O 
= ^^12/4 ^0- 

From the \\zi II7 terms, we can easily solve for £^12/4- 

^12/4 = - aj {<l>i{x), [Ho,x-Se/4\;^ <!>i{x))x 
+ a4 ($i(a;), 

- &o,2,i ( -^-sdA;!, r), <|.2(|A;|, r) )2 ( $3(y), y [i/o,, - ^10/4];' y $3(y) ) 
+ boA,o ($2(|A;|,r), $2(|A;|, r) ) 
^0 



+ 



y^2 6o,2,o//ii (2m(A;) + + 1) 



2/^2^0' 

As long as 60,4,0 7^ 0, this expression yields different values for different \k\. To see this, 
first note that the factor 



V V 2 bo,2,o AH 



2 

+ 1 ■ 



does not depend on k, and the term 

^0 



y^2 6o,2,o//ii (2m(A;) + \k\ + 1) 



2/i2l^o' 

does not depend on k. In fact, the only term that has non-trivial dependence on k in £^12/4 
is 

(^m,r), r*^m,r)}r = P + 3|fc| + P) + 6^1 + 1) ^ W + 6 .^W^ 

2 Oo,2,o A'l 

We now show that different values of k yield different values of this quantity. 

Let ki > and /c2 > be two different values of \k\ that yield the same result. Simulta- 
neously solving 

(2 + 3/ei + kl) + 6(/ci + l)m(/ei) + 6m{kif = (2 + 3A;2 + A;^) + 6{k2 + l)m(A;2) + 6m{k2f 
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and 



2m{ki) + A;i + 1 



2m{k2) + k2 + l 



forces 



m{ki) = (-3 - 5A;i + k2)/6 

m{k2) = (-3 + ki - 5k2)/Q. 

However, m{ki) and m{k2) must both be non-negative. There are no simultaneous non- 
negative solutions to 

/C2 > 3 + 5/ci 

k2 < (-3 + A;i)/5 

since this would require 3 + 5 fci < —3/5 + fci/5, which requires 24fci < —18 or ki < —3/4. 
This contradicts ki > 0, so different values of \k\ must yield different values for £^12/4. 

Therefore, at this level of perturbation, the eigenvalues generically have multiplicity 1 
when k = and multiplicity 2 when A; > 1. 



Explicitly, 




+ 




) 



ft2 

0,2,1 



{2m{k) + |A;| + 1) 



8 ^0,2,0 ^0,0,2 fJ'i 



+ 




((2 + 3|A;| + k'^) + 6{\k\ + l)m{k) + 6m{kf) 



+ 




{2m{k) + \k\ + 1). 



(3.7) 
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Order e^^/^ 

[Ho,x — £6/4] i'7/4 + [Ho,r,^ — Sg/l] "^5/4 + [Ho^y — £10/4] 1p3/4 

1 dtp 

/ /i2 -To 

+ 61,2,0 a; ^2/4 + &o,4,o ^"1/4 + 60,0,3 y^^o + ^1,0,22:2/^^0 + ^2,0,1^:^2/^0 
^0 „ , Xo / 52 9 \ 



^13/4 tpO + ^12/4 '?/'l/4- 



Order e^^/^ 

[-f^O.x — ^6/4] ^^8/4 + [Ho,r,^ — ^8/4] 'i/'6/4 + [Ho,y — S10/4] 1p4/4 

1 dlpl/4 34 2 

- 7^ rr h 03 X ^^5/4 + a4 X ?/^2/4 + 61,0,1 xy i>4/4 + 62,0,1 x yi>i/4 

I /i2 -To 'jy 

+ ^0,2,1 ^'^y V^3/4 + 61,2,0 V^3/4 + 60,4,0 r'^ 1p2/4 + 60,0,3 2/^ ^1/4 

+ 61,0,2 X ?/'i/4 + 60,2,2 ^/'o + 61,2,1 xr^y V^o + 62,2,0 V^o 

^0 „ , Xo / 92 9 \ 

—2 Ho,r,y ^2/4 + TTg ^ — — + — ^1/4 - — ^0 



+ 



2;U2Fo /i2Fo V dxdr dx J ^ ' '^^2Yi dx^ 

1 



l^2Yi 



Xo sm 7 — + — cos 7 — ^— TTT - cot e — 
or r 0*7 / \sin 0(p (77 

+ Xo cos 7 ^ 7^ 



V^o 



A4/4 V'O + ^^13/4 ^1/4 + 1^^12/4 ^2/4- 



Note: This is where we first encounter operators that mix the various different values of 
k. If we use (12 .Sp in the above expression and take ifj^ to be a linear combination of the two 
degenerate states with = A, we see that the last term on the left hand side of the equation 
contains L±i \ j, jz, A ) and L±' \ j, jz, —A ), which are linear combinations of | j, j^, A ± 1 ) 
and L±' |j, jz, —A ±1), respectively. Thus, ipQ/4 is the lowest order term that involves 
k ^ ±A. 
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Order e^^/^ 



[Hq^x — ^6/4] + [-f^0,r,7 — i'7/4, + [H^^y — £^10/4] "^b/i 

- 7, TT + «3 V^6/4 + 04 a:^ ^3/4 + 05 ^0 + &i,o,i a; y ^5/4 

+ y ^i/i + ^1,2,0 a;r^ ?/'4/4 + 6o,4,0 ^^"^ ^3/4 + Vo,3 Z/^ V^2/4 

+ &2,o,i y V^2/4 + bifi^2 xy^ tp2/4 + bo,2,2 y^ i^i/A + &i,2,i a;r^?/ ?/;i/4 

„2 „2 



+ &2,2,0 a; r + 60,4,1'^ y ^0 + ^1,4,0^'^ ^0 



-2 



+ 



/^2Vo' 



+ 



If ■ 9 Xo d 
Ao sm 7 — + — cos 7 — 



dr 



1 c> cot ^ ^ 

^7 y \ sin 6 d(j) d'-f 



X — 



+ 



Xo cos 7 sm 7 — — 

or r d'-f J od 

+ W'o 



^1/4 



+ 



d 

r sm 7 — — 
ax 



19 , (9 

777 - cot ^ — 
sm U 0(p (77 



— r cos 7 



92 



dxdO 



^15/4 V'O + ^^14/4 V'l/4 + ^^13/4 ^2/4 + ^^12/4 V'3/4- 



27 



Order e^^/"^ 

- 7; TT + «3 'ipl/A + 04 X'^ + 05 ^1/4 + ^1,0,1 Xy ?/^6/4 

z j-0 <jy 

+ ^0,2,1 ?/ ?/;5/4 + fei,2,o a;r^ ?/^5/4 + 6o,4,o '^'^ ^4/4 + &2,o,i a;^ 1/ ^3/4 

+ &0,0,3 2/^ ^/'3/4 + &1,0,2 a; ?/'3/4 + 6o,2,2 ^2/4 + &1,2,1 ^ y ^/^2/4 

+ ^2,2,0 x^r"^ V2/4 + &o,4,i r'^y ^1/4 + &i,4,o a:r^ Vi/4 + Ve.o '"^ V^o 



+ 



+ 



+ 



+ 



+ 



+ &o,o,4 I/" ^0 + bi,o,3 X y-^ ipo + b2fi,2 x"^ ipo + &3,o,i ^'^ 2/ V'o 



^0^ 



2/i2n2 
1 



Xq sin 7 



/^2ro^ 



92 



d 



dx dr dx 



^^3/4 



92 



^2/4 



9_ 

dr 



Xo d 
— cos 7 
r 



1 g cot ^ — 

d'j ) \mi 9 d(f) (97 



+ ( Xo cos 7 



d_ 
dr 



^ sin —] — 

r ^ ^ d'-f J 86 



i'2/4 



Xn 



1 



Xpy 



-f^0,r,7 ^1/4 



— r sm 7 



9_ 

9a; 



1 d 



sin 6* 



cot e 



d_ 

^7 



r cos 7 



52 



9x 96' 



^1/4 



2X02/ M ^ 92 



H2Y^ V "^a; 
1 



dx dr 



2/i2l^o' 



9 



2/^2^0' 



^0 + 



2xr — — h r — + 2x 

OT or or Ox 



/i2 V 9y 



^16/4 i'O + ^15/4 ^^1/4 + ^14/4 V'2/4 + ^^13/4 "4^3/4 + ^12/4 "^4/4- 
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3.1 The Complete Asymptotic Expansion 

We now prove the existence of a complete expansion in powers of e^^^ for the quasienergies 
and the corresponding quasimodes under suitable hypotheses. The following proposition 
completes the proof of Theorem 12.11 



Proposition 3.1 We assume the potential energy surface is smooth, with Taylor series 
given by Ii2.5\) and Ii2.6\) . Then, the eigenvalue problem for liS. 1\) can be solved by formal 
asymptotic expansions of the form 

N 

S = + 0(e(^+^)/^), 

N 

ijix, r, y, e, 0, 7) = ^ e'/^ ^i/,{x, r, y, 6, 0, 7) + 0(e(^+^)/^), 

1=0 

for any N ^N. 
Proof 

Keeping the original variables {X, R, Y), we first make use of the invariant subspace C gen- 
erated by the basis {\k)}k=-j,---j of eigenvectors of L^/, where we have dropped the fixed 
parameters j and jz from the notation. In this basis, the operator — 2L ■ J + can be 
represented by a matrix. Let I denote the identity matrix, A denote the matrix represen- 
tation of i (^—i §^ + i cos(^)^j + cos(7)^, and B denote the matrix representation of 

^ si'n(ej (""^ ^ "I" ^ cos(6')^j +sin(7)^. Note that these angular differential operators can be 
written as linear combinations of L', and L'_, which ensures that they leave C invariant. 



With these definitions, we can write 
J2 - 2L ■ J + 



dX"^ dRX dR^ \ R J dR dX 

/ V dRj R 

Then, going to the rescaled variables and dropping the symbol I, the differential operator 
3.ip takes the form 

^6/4 g2 ^8/4 ^ g2 ^ g I .2 \ e'^^' 



^ + - ^ ,Li, - 



2yUi(e) dx"^ 2/ii(e) \ Sr^ r dr r"^ J 2yU2(e) dy'^ 
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,13/4 



/i2(e) (Fo + e^/^y) dy 2 fi2{e) (Yo + e^/^y)^ [dr^ r dr 



Id 1 



+ 



/i2(e) {Yo + ey^yy 

e''/\Xo + e^/^x) 
/X2(e) (lo + e3/4y)2 

gl5/4 

/i2(e) (Yo + ey^y)^ 



92 
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dxdr dx 



or r 



,14/4 



9 , 

ox 



E 

i=2 



gl6/4 



dx'^ 



E 



r— -L^ 



1+ 



3(i+0+2fc 



x^ yK 



j + k + l>2 
k + l>l 
k even 



.1/4 



We get a matrix valued differential operator given as a formal infinite series in powers of 
by expanding the reduced masses /ij(e) and the denominators {YQ + e^^^) and {Y^ + e^/'^Y . 
Observe that in each term of the resulting expansion, the differential operators are at most 
of order two. 

The r dependence of these operators is explicit, which will allow us to check that that the 
factors 1/r and do not cause divergences in the expressions that we encounter below. 
The measure in the r variable is rdr, so the only term that might yield a vector not in 
is the L^ijr^. In the eigenspace where L^/ multiplies by zero, there is no problem. In the 
eigenspaces where L^/ multiplies by something non-zero, the wave functions contain factors 
of r, so again, there is no problem. 

We introduce the notation 



I ^i/i{x,r,y,-j) 
^i/4{x,r,y, -j + 1) 



'i/A{x,r,y) = ^ '^i/^{x,r,y,k)\k) = 

k=-j 

\ i)i/i{x,r,y,j) ) 

We have already explicitly presented perturbation theory through order e'^^ for / < 11. 
The equation we must solve at order e'/^ with / > 12 now can be expressed as 

(-f^O,x — ^6/4) ^(/-6)/4 + (-^0,r,7 — ^^8/4) ^(«-8)/4 + (-f^O,y " ^^10/4) ^(«-10)/4 

I 

+ a3X^^(Z_9)/4 + 61,0,1 3^ l/*(Z-10)/4 + ^ £'g^(/_g)/4 

g=ll 
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= Sl/i'^Qji + ■■■ + 1^^12/4 ^(i-12)/4, (3.8) 

where the symbols Dq denote at most second order differential operators in x, r, y with 
matrix valued coefficients whose entries are polynomials in these variables divided by r^, 
with p = 0, 1,2. We note also that Ho^r,j is now matrix-valued, because of the centrifugal 
term Li, /r"^, whereas iJo,a; and H^^y are scalar differential operators multiplied by the identity 
matrix. 

The point of this decomposition is to separate the vectors '^q/i of order less than or equal 
to (/ - ll)/4 from those of order (/ - 10)/4 to (/ - 6)/4. 

Let Py and Pr.,7 be the orthogonal projectors on the eigenstates $3(y) and 

on the subspace Zq = span{$2(r, \k\) \k)}k^K, respectively. We abuse notation and use 
the same symbols to denote the corresponding projectors when considered on //^(Mj., dx) ® 
L^(M+,r(ir) ® L^{9.y,dy) ® C. Note that these operators commute with one another and 
that the following identity holds for any g G N : 

a;^^+^ = P^ x'^'^^ P^, where P^ = l- Px- (3.9) 

Also, we have constructed "^1/4 so that 

^0 = PxPrnPy'^o and P^Pr^^Py-^ij^ = 0, for all / > 1. (3.10) 

Hence, for / > 1, 

*Z/4 = P^^l/^ + PxPr':^'^l/A + PxPr,-yPy^^l/4. (3.11) 

In terms of the quantities introduced in the explicit computations of the lower orders, we 
have in particular 

I 

Px'^i/, = ^ $i(x) /;/4(r,2/,fc) (3.12) 

k=-l 

PxPrn^i/i = ^ ^2{r,\k\) giii{y,k) \k) 

keK 

P^Pr,^Py^O = XI '^l(^) '^2(r,|fc|) $3(y)Cfc,o 
keK 

where Cfe^o ^ C and ^/.g^ |cfc,oP = 1- Note that by virtue of (I3.1UI] . 

9i/A{y, k) = Py giii{y, k), for any keK and any / > 0. (3.13) 
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We solve (13. 8p by two independent steps. The first consists of determining tlie vectors \l///4 
for any set of coefficients {cQ^k}k&K, and the other consists of solving an eigenvalue equation 
for Sj/4 in (Z7#(^) which may reduce the set of free coeffcients {cQ^k}k£K- It is only when we 
construct the actual quasimode that we restrict the values of the coefficients {c()^k}kGK to 
those given by the determination of the the Si/^^s. 

We now formulate our induction hypothesis for / > 12. 
IH: After solving equation (13.81) through order e^'"^-*/^ for vectors satisfying (13.101) . we have: 

• The following vectors are determined completely in terms of the coefficients {co^k}k£K 
and depend linearly on {co^k}kGK '■ 

for g = 0, 1, /-II, 

(I - P^Pr,^) ^(i_l0)/4, 

(I - ^(;_9)/4, (3.14) 

(I - P^) ^(/_8)/4, and 

(I-P.-P,^P,,,)vl/(,_7)/4. 

• The X dependence of the vector P^ Pr^^ ^(/-7)/4 is determined and has the form 

P^Pr,^^il-7y^ = P^Prn'^{l-7)/A{{9il-10)/A}), (3.15) 

with linear dependence on {g{i-io)/A{y, k)}keK, the set of functions {5'(z-io)/4} entailing 
the unknown y dependence. 

• There exist vector spaces Wq C (U*^^^ satisfying 

(Z7#W = Wo^WiD---D Wi_i (3.16) 
such that Sg/4 is determined by an eigenvalue equation in Wg, for g = 0, 1, — 1. 

Our explicit computations show that these properties are satisfied for / = 12, with Wg = 
(D*^^\ for g = 0,---,ll. We now show that the induction hypothesis holds at order e ' . 
Using (13. 9p and (I3.10p and applying P^ Pr,'yPy to equation (13.81) yields 

a3X^P^^^(i_9)/4 + bifi^ixyP;^'^(i-io)/4 + ^ P'g^'(z_g)/4 

\ 9=11 
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We note that for s = 9, 10, the vectors \l/(i_s)/4 = P^ (I — P^Pr^-y) \l/(i-s)/4 are completely 
determined by IH. By IH again, the right hand side depends linearly on the set {cQ^k}k<^K- 
Expressing the equation in the basis {$i(x)$2(|^|5 ''^)'^'3(y)}fceA' of Z2, we get a finite di- 
mensional eigenvalue equation. Restricting attention to the subspace Wi^i C C**^'^^ of free 
coefficients, we get an eigenvalue equation in Wi^i which we solve to yield Eiji and the 
subspace Wi C Wi-i of free coefficients. 

We now turn to the computation of the vectors. Application of P^ Pr^^ Py to equation 
(EH) yields 

PxPraPy'^{l-10)/A (3.17) 

= - (ifo,y - £w/A);^PxPr,^Py ia^X^P^^(^i^^y^ + 6i,o,ia;yP,^^(i_io)/4 + ^g^(Z-c)/4 

V g=ll 

where Dg = Dg — Sq/4. The right hand side is known by IH, and since Px Pr,7 -P/ ^{i-io)/4 = 
Pa; Pr,7 ^(z-io)/4) (see (13.121) . (I3.13P ). (13. lip implies that \I/(/-io)/4 is fully determined up to 
the coefficients {cQ^k}kGK- Since the dependence of Pr Pr,7 ^(«-io)/4 is linear in the previ- 
ously determined quantities, we get by IH that ^I^(/-io)/4 depends linearly in the coefficients 
{co,k}k&K- Hence, the vector P^ Pr^^ \l/(i_7)/4 ({5'(/-io)/4}) in IH is, in turn, fully determined, 
and it depends linearly on the {co^k}keK's. Thus, the same is true for (I — P^) '^{i-7)/4,- 
Application of P^ P:^^ to equation (13. 8p yields 

P.P,-;-^ ^(;_8)/4 = - (i/o,r-,7-^8/4);'PxP4 X (3.18) 

[Ho^y - ^io/4)^(/-8)/4 + a3X^P,^^(i_9)/4 + 6i,o,ia^2/^x^^a-io)/4 + Dg^i_A , 

q=ll ) 

where, by the same arguments, the right hand side is fully determined up to the coefficients 
{co,fc}fcg_ft:, on which it depends linearly. Now, from IH and the identity 

Pr^a-8)/4 = P. Pr,7 ^a-8)/4 + P.- P^ ^a-8)/4 

we see that (I — P^Pr,-^ ^(/-8)/4 is fully determined and depends linearly on the coefficients 
{co,fc}fce-ft:- 

Finally, application of P^ to equation (13.81) yields 

-Pa^^(/-6)/4 = - (-f^Cx - ^^6/4)7^ -Pa^ ^(-^0,r,7 ~ "^8/4) P^^ ^(/-8)/4+ (3.19) 

(ifo/y - ^^10/4) Pr^ ^(Z-10)/4 + a3a:^^(Z-9)/4 + &l,0,l2^^^(i-10)/4 + ^ i)g ^/--g j , 

g=ll / 
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where, this time, the right hand side is not fully determined since there is no projector 
acting on ^'(i-9)/4. However, at this step, ^'(i_io)/4 and P^^(i_8)/4 = -P^(II - PxPr^)'^ {i-fi)/A 
are fully determined and linear in the {cQ^k}k&Ki so that from IH we see that the only 
undetermined part comes from 

P,. 7 ^'y-9)/4 = ^^i{x)^2{rAk\)g(^i-si)/i{y,k)ck,o\k). 

k&K 

We conclude that the x dependence of the vector P^ ^{«-6)/4 is determined, and that the 
undetermined part of this vector depends on the set of functions {fi'(/-9)/4(?/, k)}k(^K purely 
linearly. 

Thus, we have reproduced the all the requirements of the induction hypothesis, which 
ends the proof. H 

3.2 The Expansion Around a Local Minimum 

We now describe the construction of quasimodes of arbitrarily high order under assumptions 
that are only local. This construction uses the formal expansions of Proposition 13. II and the 
insertion of cutoff functions. The construction is quite similar to that given in [5], so we 
refrain from presenting all details. 
Let > be fixed and set 

N 

vi/(^)(x,r,y,^,0,7) = J] e'/4^,/4(a;,r,y,^,0,7), 

1=0 

N 

^iN) ^ ^ ^//4 ^^^^^ (3 20) 

V^''\X,Y,R) = Yl MX-XoY + e J] b.^kA^ - X^Y {Y - YoY , 

l<{N+l)/3 j + k + l>2 

k + l>l 
k even 
4 + 3{j + l) + 2k< N 

where the vectors 'i/'z/4 and the scalars Si/ 4^ are defined in Proposition 13. 1[ 

Then we introduce a cutoff function. Let T : M ^ [0,1] be C°° and such that supp c 
[-2,2] with J^{t) = 1 for t e [-1,1]. We set 

^,(X,P,F) = J^{{X - Xo)/e'^) J^{R/e'') J^{{Y -Yo)/e''), 

where < 61 < 3/4, < ^2 < 1/2 and < ^3 < 3/4. 
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The quasimode is defined as 

vi>g')(x,i?,y,^,0,7) 



(3.21) 



= e-^/^ Tlx, R, Y) vl/W((X - Xo)/e'/\ R/e''\ {Y - Y,)/e''\ 6, 7). 

The factor of e~^/^ in this expression ensures asymptotic normahzation of the quasimode 
because of the Jacobian factor in the integral for the norm. 



Proposition 3.2 Let 



Hie) 



Vi(X) + e V^iX, R, Y), 



satisfy the hypotheses of Proposition [X71 Then, for any G N, there exists a constant C^, 
such that the vector ^3.21\) and the scalar ^3.20^ satisfy \\^^q^\\ = 1 +0(e^/^) and 



(TV) 



(iV) 

Q 



< Cm e(^+^)/^ as e ^ 0. 



Proof 

We begin by computing the norm of "^^q^ ■ The vectors ipi/^,., for / = 0, ■ ■ ■ , A^, are given as 
a finite hnear combinations of angular functions \k,jz,j), [k = —j,---,j), multiplied by 
Gaussians in x,r,y, times polynomials in these variables. Thus, they all belong to L^. 
In particular, by our choices for tpQ, we have 

J |e-^/^ Mix - Xo)/e3/^ R/e^/\ (Y - Yo)/e^/\ 6, 0, 7) ^ RdRdX dY dn 



\ipo{x,r,y,9,(l),'y)\'^ rdrdxdydVt 



where dVL denotes the solid angle element in the angular variables. The norms of the other 
tpiH are similarly 0(1). 



Hence 



W||2 
Q II 



|^W + (jr2_i)^(iV)||2^ where. 



'1_^2^^(7V)| 



(3.22) 



< 



|^(^)((X - Xo)/e3/^ i^/e^/^ {Y - Y^)/e'"\ 9, 4>, ^)\^ RdRdX dY dQ. 



R > e^^ 
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The choice of exponents 5j and the exponential decay of \1>(^) imply that fl3.22p is of order 
e°°, and we finally see that 

= 1 + 0(ei/^). 



By construction, there exist C > and -D > 0, independent of e, such that 
7^(^)(X,i^,r) = Vi{X) + eV2{X,R,Y) - V^''\X,Y,R) 

satisfies 

\n'^^\X,R,Y)\ < C(|X-Xo|(^+^)/3 + e|X-Xo|'^i?''|F-Fo|'), (3.23) 
where 4 + 3(a + c) + 2c > + 1, if (|X - Xo| + i? + |r - Yq\) < D. Consider now 

Due to the support conditions imposed by the cutoff, we can estimate J-'^TZ^'^^ by means of 
fl3.23p . and, after passing to the rescaled variables x, r, y, we obtain 

J^,{X,R,Y)\n^^\X,R,Y)\ < J-,(X,/?,F)e(^+^)/'^C {\x\^^+^^/'' + \x\''r''\y\') . 

Once again using the Gaussian decay of we finally get the estimate 



We now have estimated everything except the terms in which the kinetic energy acts on 
the cutoffs. First note that derivatives with respect to angular variables do not affect the 
cutoffs. Next, by the Leibniz formula, the first and second derivatives with respect to x, y, 
or r acting on JF^ v[/(^) yield supplementary terms given by first and second derivatives of JF^ 
multiplied by vl>(^) or first derivatives of \1>(^). By construction of the cutoff, the successive 
derivatives of JF^ are supported away of the origin in at least one of the variables x, y, or 
r. Since and its derivatives are Gaussian times polynomials in these variables, these 
supplementary terms are all of order e°°. 

Finally, taking into account the formal expansions of Theorem 13.11 and the definition 



2/ii(e) 



A 



2/i2(e) 



A(n,y„y3) + V^^'^X, R, Y,), 
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we get the norm estimate 



0(e°°) 



4 Inclusion of the Electrons 

In this section we show that including the quantum mechanical treatment of the electrons 
does not change the expression for the energy up to an error of order e'^. 

We decompose the Hamiltonian for all the particles in the molecule as the sum of the 
nuclear kinetic energy plus a self-adjoint electron Hamiltonian hi{Y, 6, (p, R, 7, X). The 
electron Hamiltonian depends parametrically on {Y, 6, 0, R, 7, X) and acts on functions of 
all of the electron variables, that we describe jointly with the single symbol Z. To avoid 
questions about Berry phases, we assume hi(Y, 9, 0, R, 7, X) commutes with complex con- 
jugation, i.e., it is a real symmetric operator. 

Because of rotational symmetries, the electron Hamiltonian can be written as 

h^{Y, 9, 0, R, 7, X) = U{9, 0, 7) hiX, R, Y) U{9, 0, 

where U{9, 0, 7) is unitary on the electron Hilbert space and depends smoothly on 9, </>, and 7. 
As a consequence, discrete eigenvalues of hi{Y, 9, 0, R, 7, X) do not depend on 9, 0, or 7. 

We assume that the resolvent of h2{X, R, Y) depends smoothly on {X, R, Y). As a 
result, all discrete eigenvalues of hi(Y, 9, 0, R, 7, X) depend smoothly on the nuclear con- 
figurations. 

We assume further that the ground state eigenvalue V{X, R, Y) of h(Y, 9, 0, R, 7, X) is 
discrete and non-degenerate for each fixed value of [Y, 9, 0, R, 7, X). We also assume that 
V{X, R, Y) has a global minimum at (Xq, 0, ^0) with a strictly positive Hessian at that 
minimum. To ensure that we are approximating discrete eigenvalues for the full molecular 
Hamiltonian, we assume that the V{Xo, 0, Yq) is strictly below the bottom of the spectrum 
of /i2(X, R, Y) for all (X, R, Y) outside a small neighborhood of (Xq, 0, Yq). 
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We now introduce e-dependence in h2, and hence hi. We choose functions Vi{X) and 
V2{X, R, Y) that satisfy 



V{X, R, Y) = Vi{X) + eo V^iX, R, Y) 

and the restrictions imposed after expression (12 .4^ . Here eo is a fixed value of e that we 
take to be the fourth root of the electron mass divided by the carbon C^^ nuclear mass. We 
then define h{e, Y, 6, 0, R, 7, X) by replacing V{X, R, Y) by Vi{X) + e V2{X, R, Y) in 
the spectral decomposition of hi{Y, 9, 0, R, 7, X). Thus, we only introduce e-dependence 
in this single eigenvalue and alter none of the eigenfunctions. 

Remark To minimize technicalities, we have made assumptions for all {X, R, Y). At the 
expense of inserting cut off functions, our assumptions need only be imposed for [X, R, Y) 
in a neighborhood of {Xq, 0, Yq). 

We shall write down an explicit quasimode with an 0(e^^/^) energy error for the Schro- 
dinger operator 

e^ e^ 
-^(^) ^ ~ o — r\ ^(^1,^2,^3) - 7^ — rr- ^{¥1X2X3) + ^(^' Xi, X2, X3, Yi, Y2, Fs), 

rewritten in terms of the variables [Y, 6, </>, R, 7, X, Z). 
The quasienergy will be 

S{e) = So + e'l' £,i, + e'" + e'"' S^o/a. (4.1) 

but the quasimode will be somewhat complicated. 

To specify the quasimode, we first let x{^i 'Pi Ri I1 X, Z) denote a normalized real 
ground state eigenfunction of /i(e, y, 9, (p, R, 7, X) that depends continuously on its vari- 
ables. Next, we let 



C(e, Y, 9, 0, R, 7, X) = e-/^ fl ^'/4 (^^, J,, 0, 0, 7 



where the ipi/i the wave functions from Section [3] with g2/4{y, ±A) = gz/iiy, ±A) = 
fi/i{r, y, k) = /5/4(r, y, k) = 0. Note that when A = there is one hnearly independent 
choice for (. When A > 0, we have two linearly independent choices corresponding to 
k = ±A. 
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The quasimode is 
vl/(e, F, e, 0, R, 7, X, Z) 



i?, y) C(e, 0, i?, 7, X) xiY, Q, 0, i?, 7, X, Z) 



+ -- Tlx, R, Y) h{e, Y, 6, 0, R, 7, X) - V{e, X, R, Y) 



-1 



(4.2) 



X I ^(6, Y, e, 0, i?, 7, ^) ^, 0, R, 7, ^, ^) 

+ ||(e, F, ^, 0, 7, X) ^, 0, R, 7, X, Z) 



Theorem 4.1 There exists a constant C, such that the function \l/(e) given by ^4-2^ and 
quasienergy S{e) given by ^.1] ) satisfy 

II ^(e) II = 1 + O 

and 



Hie] - £{e] ^fe, ■ 



(4.3) 



Proof The function \I/(e, ■ ) equals the normahzed vector ipQ x plus terms that are 
orthogonal to -00 X- Since the largest of these orthogonal terms is e^^^ -01/4 x, we see 
that \I^(e) has norm 1 + 0(6^/^). 

To prove the second estimate of the theorem, we begin by noting that the electronic 
eigenfunction x has the form 

X(F, <\>, R, 7, X, Z) = U{9, 0, 7) Xo{Y, R, X, Z), 

where U{9, 0, 7) is unitary. 
We next compute 

(i7(6) - ^(e)) T{X, R, Y) C(e, ■) x(-), (4-4) 
where H{e) is decomposed as 



Hie] 



A 



(X-i,X2,Xz) 



2^i(e) — 2/i2(e) 



+ 



/i(e, Xi, X2, X3, Fi, 1^2, Y^) - V^{X) - e V^iX, R, Y) 



+ Vi{X) + e \/2(X, i?, F), 
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with the two final terms expanded in their Taylor series of appropriate orders. We write the 
resulting expression in the variables {Y, 6, 0, R, 7, X, Z). When the all the derivatives in 
H{e) act on (, all terms that are larger than order cancel because of Taylor series estimates 
and the choices of the When all the derivatives act on x, all terms are 0{e^) or smaller 
because x is smooth and the cutoffs are zero the singularity at F = 0. When any derivatives 
act on JF^, we obtain terms of order 0{e'^), for any q, due to the rapid fall off of the functions 
in (. The term that arises from [h{e) — Vi — e V2] yields zero because it acts only on the x- 
The remaining terms in (14.41) contain terms in which a partial derivative acts on ( and 
the same partial derivative acts on x- All of these terms are O(e^) or smaller, except for 

^(e, Y, 9, 0, R, 7, X) 9, 0, R, 7, X, Z) (4.5) 

+ ||(e, Y, 9, 0, R, 7, X) ||(F, 9, 0, R, 7, X, Z). 

Thus, (14. 4p yields (14.51) plus O(e^). However, when the [h{e) — Vi — €¥2] acts on the second 
term in (14.21) . the terms that arise from (14.51) cancel, leaving us with 0{e^) errors plus the 
kinetic energy and potential terms acting on the second term in (14. ip . Because of the cutoff, 
the potential terms yield bounded operators times O(e^) terms. When the kinetic energy 
acts on these terms, we obtain terms of order e^/^ or smaller, since everything is smooth, 
and the largest terms come from and two X-derivatives acting on (. 

dx 

Note that when computing the norm in (14. 3p . it is essential that x be orthogonal to —— 

oX 

and 7- — , or cross terms would yield terms of order greater than e^. This orthogonality is 
oR 

guaranteed by our hypothesis that the electron Hamiltonian /i(e, ■ ) be real symmetric and 
that we choose x to be real. H 
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